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ABSTRACT 

We investigate the possibility of habitable Trojan planets in the HD 23079 
star-planet system. This system consists of a solar-type star and a Jupiter-type 
planet, which orbits the star near the outer edge of the stellar habitable zone in an 
orbit of low eccentricity. We find that in agreement with previous studies Earth- 
mass habitable Trojan planets are possible in this system, although the success 
of staying within the zone of habitability is significantly affected by the orbital 
parameters of the giant planet and by the initial condition of the theoretical 
Earth- mass planet. In one of our simulations, the Earth-mass planet is captured 
by the giant planet and thus becomes a habitable moon. 

Subject headings: extrasolar planets, habitable zone, orbital stability and plan- 
etary climate. 

1. Introduction 

The existence of planets in orbit about solar-type stars is now a well-established ob- 
servational result. Obviously, the ultimate quest of these studies is to discover Earth-hke 
planets located in the habitable zones (HZs) of their host sta rs. So far, a smal l number of 



super -Earth planets with masses of up to about 12 (e.g., lUdry et a/.l 120071 : IVogt et al. 



20101 ) have been found, typically located around M-type dwarf stars. Nevertheless, the exis- 
tence of Earth-mass planets, including those hosted by solar-type stars, is strongly implied 
by various observational findings inclu ding the occurrence and mass distribution of close-in 



super-Earths, Neptunes, and Jupiters (IHoward et a/.ll2010l ). Measurements by the authors 
indicate an increasing planet occurrence with decreasing planetary mass Mp akin to M~^'^^, 
implying that 23% of stars harbor a close-in Earth-mass planet (ranging from 0.5 to 2.0 M^); 
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see also iMarcy fc Butlerl ( 120001 ) for earlier results. Very recent support for the existence of 
Earth-type planets outside the Solar System is lent by the disc overy of Kepler- 10b, Kepler's 
first rocky planet with an estimated mass of 4.6 Earth masses (IBatalha et a/.ll201ll ). 



Long-term orbital stability of Earth-like planets in stellar HZs is necessary for the evo- 
lution of any form of life, particularly intelligent life. There is a large array of studies 
focusing on the orbital stability of hypothetical Earth-mass planets in stellar HZs concern- 
in g different type s of host stars and star-p l anet config u rations. Examples in c lude studies 



Gehman et al 




Jones et al. 



Menou fc Tabachnik J2003h 



(12001 



2005 



Cuntz et al 



20061) IJones fc Sleepl (120101) . iNoble et al 



J2003h . Ivon Bloh et all N2003h . and 



Asghari et al 



Particular types of systems are those where a Jupiter-type planet orbits a star in 
the stellar HZ, therefore jeopardizing the possibility of habitable terrestrial planets in that 
system. This is actually the situation of HD 23079, the focus of the present paper. 



Previously, iNoble et all (120021 ) investigated the orbital stabihty of terrestrial planets 
inside the HZs of 47 UMa and HD 210277. The center stars of these systems are very similar 
to the Sun concerning mass, spectral type, and effective temperature. Orbital stability was 
attained for the inner part of the HZ of 47 UMa; however, no orbital stability was found 
for hypothetical Earth-mass planets in the HZ of HD 210277. In this case, a Jupiter-type 
planet crosses the stellar HZ, thus effectively thwartin g habitability for this system. Very 
recent examples were also given by lYeager et all (120111 ) who studied the star-planet systems 
HD 20782 and HD 188015. In both cases, the giant planet significantly interferes with any 
Earth-mass planet in the stellar HZ assumed to have formed there for the sake of study. In 
all cases, the Earth-mass planet was ejected from the stellar HZ in a very short time. 

However, if a giant planet is orbiting the star in the stellar HZ, there is stil l the principal 
possib ility of habitable Tro j an pl anets in those systems as pointed out by, e.g., [Dvorak et al 
(I2OO4J ) and ISchwarz et all (120071 ). A Trojan planet is one located around one of the La- 
grangian points L4 and L5 of the giant planet. These points lie on the giant planets 
orbit, ahead (L4) and behind (L5) the planet, each forming an equilateral triangle with 
the plan et and its sta r. Thu s, Trojan planets are also in a 1:1 resonance with the giant 
planet. [Dvorak et al\ (120041 ) investigated the stability regions of hypothetical terrestrial 
planets around L4 and L5 in specific systems, including HD 23079, in the framework of the 
restricted three body problem. They obtained relationships between the size of the stability 
regions and the orbital parameters of the gi ant planets, p a rticul arly its eccentricity. Studies 
about Neptune Trojans were given by, e.g., [Dvorak et all (120071 ). 



A study by ISchwarz et al\ (120071 ) identified several exoplanetary systems that can har- 
bour Trojan planets with stable orbits in the stellar HZs. Concerning HD 23079, this study 
concluded that a Trojan planet will only spend 35% of its time in the stellar HZ, assumed 
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to extend from 0.85 to 1.60 AU. In our st udy, however, we y y ill co nsider a zone of habit- 
ability based on the generahzed estimate by lUnderwood et al\ (120031 ) . implying habitability 
between 0.99 and 1.97 AU (see below). This means that the habitable area, w hich is the area 
of the stellar HZ annulus, is increased by 58% compared to that considered by lSchwarz et al. 
teOOTl ). In the present study, we will conclude that habitable Trojan planets are indeed 
possible in the system of HD 23079, although their existence is significantly affected by, e.g., 
the orbital parameters of the giant planet. Next we will describe our theoretical approach. 
We will discuss the adopted methods and the system parameters for HD 23079. Thereafter, 
we will describe our results. Finally, we will present our conclusions. 



2. Theoretical Approach 
2.1. Stellar and Planetary Parameters 



HD 23079 has been mo nitored as part of the Anglo-Australian Planet Search (AAPS) 
program (ITinney et a/.ll2002l ) that is able to perform extrasolar planet detection and measure- 
ments with a long-term, systematic radial velocity precision of 3 m s~^ or better. HD 23079 
was identified to host a Jupiter-type planet in a relatively l arge and nearly circular or- 
bit. HD 23079 is an inactive main-sequence star; iGray et al\ ( 120061 ) classified it as F9.5 V 
(see Table 1; al l para meters have their usual meaning), an updated result compared to 
Houk fc Cowleyl (Il975l ) who found that HD 23079 is intermediate between an F8 and GO 
star. Its stellar spectral type corresponds to a mass of M = 1.10 ± 0.15 Mq. The stellar 
effective tem perature and radi us are given as Teg = 6030 ± 52 K and R = 1.106 ± 0.022 Rq, 
respectively ( jRibas et a/.ll2003l ). Thus, HD 23079 is fairly similar to the Sun, though slightly 
hotter and slightly more massive. The detected planet (HD 23079 b) has a minimum mass 
of Mpsini = 2.45 ±0.21 Mj. Furthermore, it has a semimaior axis of = 1.596 ±0.093 AU 
and an eccentricity of Cp = 0.102 ± 0.031 ( Butler et al. 2006h. correspondin g to an orbital 
period of P = 730.6 ± 5.7 days. The original results by lTinney et al\ (|2002[ ) indicated very 
similar planetary parameters. 

The orbital parameters of HD 23079 b are relatively similar to those of Mars, implying 
that HD 23079 b is orbiting its host star in or near the outskirts of the stellar HZ; see 
discussion below. The existence of HD 23079 b, a planet even more massive than Jupiter, 
makes it difficult for a terrestrial planet to orbit HD 23079 at a similar distance w ithout being 



heavil y affe cted by the g i ant pl anet; see results from previous case studies by iNoble et al. 



(]2002[ ) and lYeaeer et al\ (1201 If ) who focused on the dynamics of HD 20782, HD 188015, 
and HD 210277. Concerning HD 230 79, a previous inves tigation pertaining to habitable 
terrestrial Trojan planets was given by [Dvorak et all ( 120041 ). 
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2.2. Method of Integration 



For our simulations of the HD 23079 system, we consider both the observed giant planet 
and a hypothetical terrestrial planet of one Earth-mass, i.e., 3.005 x 10~^ M©, which allows 
us to execute a grid of model si mulations. T he method of integration uses a fourth-order 
Runge-Kutta integration scheme (lGarciall2000l ). The code has been extensively tested against 
known analytical solutions, including the two-body and restricted three-body problem (see 



Noble et a/.ll2002l : ICuntz et a/.l 120071 : lEberle et a/.ll2008l . for detailed results). In the frame- 
work of our simulations that we limit to 10^ yrs, we apply a time-step of 10~^ yrs for the 
integration scheme that is found to be fully appropriate. In that regard, we pursued test 
studies comparing the planetary orbits based on three different integration time-steps, which 
are: 10~^, lO"'' and 10~^ yrs. In particular, we evaluated ARij, i.e., the magnitude of the 
difference between the position of the planet when different step sizes of 10~* and 10"-' 
were used. We found that there is no significant change in outcome between models with 
time-steps of 10~^ and 10"^ yrs. 

The initial conditions (i.e., starting velocities) for the orbits of the Earth- mass planets 
were chosen such that the planet was assumed to start at the midpoint of the stellar HZ 
(1.4779 AU) and to be in a circular orbit about the star, although it is evident that it will be 
significantly affected immediately by gravitational pull of the giant planet, which will prevent 
the planet from continuing a circular motion. For each set of models, defined by sets of values 
for the semimajor axis and eccentricity Cp, given as 1.503, 1.596, 1.689 AU and 0.071, 
0.102, 0.133, respectively, 8 different configurations are considered. They are defined by the 
8 different starting (phase) angles for the Earth- mass planet, which are varied in increments 
of 45° noting that 0° corresponds to the 3 o'clock position. Moreover, the starting position 
of the Jupiter- type planet (HD 23079 b), for which we assume its minimum mass value of 
2.45 Mj, was varied between its periastron and its apastron position. Therefore, a total of 
144 initial configurations has been considered. Note that the Jupiter-type planet was always 
started at the 3 o'clock position, which after adjusting the orbital layout of the giant planet 
always coincided with its periastron (see Fig. 1) or apastron position depending on the type 
of model. 



2.3. Stellar Habitable Zone 



The extent of the HZ of HD 23079 has been calculated fo l lowing the formalism by 



Underwood et al\ (120031 ) based on previous work by iKasting et al\ (119931 ) . lUnderwood et al. 
( 120031 ) supplied a polynomial fit depending on the stellar luminosity and the stellar effective 
temperature that allows to calculate the extent of the conservative and the generalized HZ. 
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Noting that HD 23079 is more luminous than the Sun, it is expected that its HZ is more 
extended than the solar HZ, for which t he inner and outer l imit of the generalized HZ were 
given as 0.84 and 1.67 AU, respectively (IKasting et a/.lll993l ). The generalized HZ is defined 
as bordered by the runaway greenhouse effect (inner limit), where water vapour enhances the 
greenhous effect thus leading to runaway surface warming, and by the maximum greenhouse 
effect (outer limit), where a surface temperature of 273 K can still be maintained by a cloud- 
free CO2 atmosphere. The inner limit of the conservative HZ is defined by the onset of water 
loss, i.e., the atmosphere is warm enough to allow for a wet stratosphere from where water 
is gradually lost by photodissociation and subsequent hydrogen loss to space. Furthermore, 
the outer limit of the conservative HZ is defined by the first CO2 condensation attained by 
the onset of formation of CO2 clouds at a temperature of 273 K. 

For HD 23079, the limits of the conservative HZ are given as 1.1378 and 1.6362 AU, 
whereas the limits of generalized HZ are given as 0.9896 and 1.9662 AU (see Fig. 1). The 
limits of the generalized HZ are those employed in our numerical planetary studie^j. The 
underlying definition of habitability is based on the assumption that liquid surface water is 
a prerequisite for life, a key concept that is also the basis of ongoing and futur e searches 
for extrasolar habitable planets (e.g., ICatanzarite et a/.ll2006l : ICockell et a/.l 120091 ). The nu- 
merical evaluation of these limits is based on an Earth-type planet with a CO2/H2O/N2 
atmosphere. Specifically, the inner limit of habitability is set by the loss of water from the 
upper planetary atmosphere through photodissociation and subsequent escape of hydrogen 
to space associated with a run- away greer ihouse effect. The outer limit of habitabili ty is 
given by the maximum greenhouse effect (IKasting et a/.lll993l : lUnderwood et a/.l 120031 ) . by 
which a surface temperature of 273 K can be maintained by a cloud-free CO2 atmosphere. 

We point out that concerning the outer e dge of habitability, even less conservative limits 



have been proposed in the meantime (e.g., iForget fc PierrehumbertI 119971 : iMischna et al. 



2OOOI ). They are based on the assumption of relatively thick planetary CO2 atmospheres and 
invoke strong backwarming that may further be enhanced by the presence of CO2 crystals 
and clouds. However, as these limits, which can be as large as 2.4 AU in case of the Sun, 
depend on distinct properties of the planetary atmosphere, they are not relevant for our 
study. Nevertheless, we convey this type of limit for the sake of curiosity (see Fig. 1), noting 
that it has properly been adjusted to 2.75 AU in consideration of the radiative conditions 
of the planetary host star, HD 23079. Moreover, the significance of th is extreme lirnit ha s 
recently been challenged based on detailed radiative transfer simulations (IHalevy et a/.ll2009l ) . 



^The physical limits of habitability are much less stringent than implied by the numerical precision of 
these values; nevertheless, these values were used for checking if the Earth-mass planet has left the stellar 
HZ. 



- 6 - 



3. Results and Discussion 

3.1. Case Studies of Habitable Trojan Planets 

Table 2 and Table 3 summarize the time the Earth-mass planet remains within the 
stellar HZ, i.e., before exiting the stellar HZ or being permanently ejected from the system. 
Of the 144 total considered initial configurations 13 survived at least 1 million years, 93 
crossed the upper limit of the HZ, 28 crossed the lower limit of the HZ, and 10 collided or 
may have had a very close approach with the giant planet. Some of those who crossed the 
HZ at the lower or upper limit as first exit from the HZ may have had a very close approach 



with the giant planet, resulting in destruction while entering the Roche limit (jWilliamsll2003l ) 
or in a collision with the giant planet, at a later time. Of the 13 survivors, 12 are Trojan 
types, that is they exist i n stable orbits aroun d the equilateral equilibrium positions much 



like that demonstrated in [Dvorak et al\ (120041 ) . 



In the cases where the giant planet is initially in the periastron position, only models 
with the smallest considered semimajor axis and eccentricity combination, which are ap = 
1.503 AU and Cp = 0.071, result in habitable Trojan planets (see Fig. 2). In this case, four 
different starting positions (phase angles) appear to be consistent with long-term stability 
(see Table 2). It is clear that the Earth- mass planet is safely inside of the stellar HZ but it 
is a snug fit. For the next larger eccentricity considered, which is 0.102, there are various 
cases where the Earth-mass planet stays within the HZ for some hundreds thousand years 
before finally crossing the upper limit of the HZ. When the eccentricity of the giant planet 
is increased to 0.133, the Earth-mass planet remains within the HZ at best for only a few 
hundred years. 

The situation is, in principle, similar for the cases where the giant planet is initially 
placed at the apastron position. In this case, for ap = 1.503 AU, the Earth- mass planet 
remained in the stellar HZ for at least a million years for two eccentricity simulations, which 
are Cp = 0.071 and 0.102 (see Figs. 3 and 4, respectively). Comparing Fig. 4 to Fig. 3, it is 
clear that the Earth-mass planet moves in a wider area and approaches the edges of the HZ 
for the larger eccentricity, thus illustrating how the planet remained in the HZ for such a 
long time before exiting in the periastron case with the same parameters. In some of those 
latter cases, we found that the planet was outside the HZ for a brief time (i.e., considerably 
less than a planetary orbit), but most li kely without losing i t s hab itability. This conclusion 



is motivated by the previous study of IWilliams &: PoUardl ( 120021 ) that showed that brief 



excursions from the HZ are insufficient to nullify planetary habitability because the latter is 
expected to mainly depend on the average stellar ffux received over an entire orbit, rather 
than the length of the time spent within the HZ. 
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Note that Figs. 2 and 3 display models of orbital stability for the Earth-mass planet in 
a synodic (rotating) coordinate system. Thus, the "banana-shaped" areas correspond to the 
domains about L4 or L5, where stability for the Earth- mass planet is encountered. The thin 
line at the 3 o'clock position corresponds to the motion of the giant planet due to its slightly 
elliptical orbit. Clearly, only Earth- mass planets placed at phase angles of 45°, 90°, 270°, 
and 315° have a reasonable chance to develop into Trojan planets, whereas for other starting 
angles ejections from the HZ, and usually also from the star-planet system, will occur due to 
gravitational interaction with the giant planet. If the Earth-mass planet was initially placed 
at an angle of 60° or 300°, it can be expected that it will continue to remain a Trojan planet. 

For the sake of curiosity, we also evaluated various cases where the Earth-mass planet 
never had a chance of becoming habitable. Hence, we chose five cases of different semimajor 
axes and eccentricities for the giant planet. In all cases the giant planet started at the 
periastron position and the initial phase angle of the Earth-mass planet was chosen as 180°. 
The simulations are depicted in Fig. 5. In Fig. 5a, with ap = 1.503 AU and Cp = 0.102, the 
system experiences a relatively long period during which the Earth-mass planet first exits 
the HZ at 48.4 yrs. This event is preceded by a close approach with the giant planet. The 
simulation is terminated at 195.9 yrs due to an expected collision with the giant planet. 

Figures 5b to 5d are all based on ap = 1.596 AU, but the depicted simulations assume 
different eccentricities for the giant planet, which are Cp = 0.071, 0.102, and 0.133, respec- 
tively. In Fig. 5b, the system experiences a short period during which the Earth-mass planet 
first exits the HZ at 7.93 yrs. This event is again preceded by a close approach with the 
giant planet. The simulation is terminated at 9112 yrs due to the expected collision with the 
giant planet. In Fig. 5c, the system experiences a short period during which the Earth- mass 
planet first exits the HZ at 7.80 yrs. This event is preceded by a close approach with the 
giant planet. The simulation is terminated at 13400 yrs considering that the Earth-mass 
planet is ejected from the system. Habitability is ultimately prevented as the Earth-mass 
planet becomes "free-floa ting". Free-floating p lanets have previously been observed in case 
of the Trapezium cluster (jhucas fc Roche 2000 ): note that planetary ejections due to orbital 



instabilities are an important candidate process for this finding. 

In Fig. 5d, with ap = 1.596 AU and Cp = 0.133, this system experiences a short period 
during which the Earth-mass planet first exits the HZ at 7.82 yrs. It is reentering and exiting 
the HZ several times. However, the simulation is terminated at 54.60 yrs due to an expected 
collision with the giant planet. In case of Fig. 5e, with ap = 1.689 AU and Cp = 0.102, the 
system again experiences a short period during which the Earth-mass planet first exits the HZ 
at 4.78 yrs. This event is preceded by a close approach with the giant planet. Eventually, the 
planet also becomes free-floating; the simulation is terminated at about 7 x 10^ yrs. Figure 
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5a and 5d show oscillatory behaviours regarding the orbital motion of the Earth-mass planet. 
Noting that the Earth-mass planet starts at a phase angle of 180°, it initially orbits the star. 
However, when it approaches the giant planet, its orbit is being perturbed causing the loops. 
Thus, the Earth-mass planet exits and re-enters the HZ multiple times until the end of the 
simulation. 



3.2. On the Possibility of Habitable Moons 



Our set of model simulations reveal a considerable variety in the dynamics of the Earth- 
mass planet. The most surprising case is the following: For = 1.596 AU and Cp = 0.133 
(see Fig. 6) with the Jupiter-type planet initially placed at periastron position and the Earth- 
mass planet placed at 0°, it was found that the latter never crosses the inner or outer limit 
of the stellar HZ during the simulation time of 10^ years. However, it is found to orbit the 
giant planet in a retrograde orbit (relative to the orbital motion of the giant planet about 
the star). In this case, the Earth- mass planet is captured by the giant planet and becomes 
a habitable moon, which occurs almost immediately after the start of the simulation. 

The analysis of its orbital data shows that the moon's semimajor axis concerning its 
motion about the giant planet is amoon — 0.051 AU. Its eccentricity is Cmoon — 0.8 entailing a 
perigee and apogee of 0.0034 and 0.098 AU, respectively. Thus, with a uniform data sampling 
rate, the moon is most likely to be recorded at or near apogee. From Fig. 7 it is evident that 
there is also a precession of the perigee in a retrograde sense with a period of approximately 
30 years. Figure 8 shows two histograms regarding the time-dependent distance of the moon 
from the giant planet, which reconfirms the moon's highly eccentric orbit. The existence of a 
habitable moon in the HP 23079 system is also consistent with the criterion of Hill stability 
as pointed out by, e.g., iDonnisonI (120 lOl ). This study explores dynamic Hill stability for a 
large variety of three-body systems considering moon/planet mass ratios of 0.1, 0.01 and 
0.001. 



There is a persistent interest in the study of habitable moons with respect to extrasolar 
giant planets orbiting host stars in the stellar H Zs. Previous s t udies o f habitable moons 
i n sys tems akin to HD 23079 have been given bv IWilliams et all (jl997l ). iBarnes fc O'Brien 



(I2OO2I), and others. The study by IWilliams et al\ (119971 ) did not include HD 23079b as 
this star-planet system was unknown at the time when this study was pursu ed. However, 
by targeting the companions of 16 Cyg B and 47 UMa, IWilliams et al\ ( 119971 ) investigated 
appropriate orbital parameters of possible moons, and pointed out that the moons need to 
be large enough (i.e., > 0.12 M®) to retain a substantial and long-lived atmosphere, and 
furthermore would need to possess a significant magnetic field to prevent its atmosphere 
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from being sputtered away by the ongoing bombardment of energetic ions from the planet's 
magnetosphere, if existing. Another study of pos sible moons, which i s fully applicable to the 
HD 23079 star-planet system, has been given by lBarnes fc O'BrienI (120021 ) . They concluded 
that Earth-like moons of a Jovian planet like HD 23079b would be able to exist for at least 
5 Gyr considering that the stellar mass of HD 23079 exceeds 0.15 Mq. 



4. Conclusions 

The aim of our study was to add to the investigation of habitable Trojan planets in 
the HD 23079 star-planet system. This system consists of a main-sequence star slightly 
hotter than the Sun. Additionally, it contains a Jupiter-type planet with a minimum mass 
of 2.45 Mj that is orbiting the star in a slightly elliptical orbit that is positioned within 
the stellar HZ. The main goal of our study was to explore if Earth-mass habitable Trojan 
planets can exist in this system. 

As the centerpiece of our study, we calculated a total of 144 orbital stability simula- 
tions for the Earth-mass planet by choosing different starting positions (phase angle) as well 
as placing the Jupiter-type planet either at periastron or apastron position. The attain- 
ment of habitability solutions was found to critically depend on various parameters, which 
include the orbital parameters of the giant planet (semi-major axis, eccentricity) and the 
initial condition (phase angle) of the theoretical Earth- mass planet. We encoutered a va- 
riety of different outcomes, which include (1) ejection of the Earth-mass planet from the 
system, (2) engulfment of the planet by the star (or possible destruction in accord with 
the Roche limit criterion), (3) capture of the planet, thus becoming an habitable moon, or 
(4) remaining within the stellar HZ. The latter case was only attained in models where the 
orbit of the giant planet had a relatively low eccentricity (but still within its observation- 
ally given uncertainty), which however may be partially due to the implemented choice of 
planetary starting positions. Concerning the latter case, there were also cases (not shown 
in detail) where the planet took short-term excursions from the HZ (i.e., considerably less 
than the orbital period of HD 23079b, which is about 730 d), which should be insufficient 
to nullify its habitability because the latter is expected to mainly depend on the average 
stellar fl ux received over an entir e orbit, rather than the length of the time spent within 



the HZ (jWilliams fc PollardI 120021 ). although the ultimate effect of temporarily leaving the 



zone of habitabili ty will still partially depend on the atmospheric thickness, structure and 



composition (e.g.. [Dressing et a/.ll2010l ) 



Moreover, we note that our study is supplementing previous work by ISchwarz et al 



( 120071 ) who concluded that a Trojan planet in the HD 23079 star-planet system will only 
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spend 35% of its time in the stellar HZ. However, this estimation was based on a considerably 
narrower zo ne of habitabili t y tha n used in the present study. Another, albeit minor, differ- 
ence is that ISchwarz et all (120071 ) used slightly different orbital parameters for HD 23079b 
than in the current study. In conclusion, it can be argued that the system of HD 23079 is 
very well suited for the existence of habitable Earth-type Trojan planets, and thus deserves 
serious consideration in ongoing and future planetary search missions. 
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Fig. 1. — Extent of the HZ for HD 23079, defined by its conservative limits (dark grey) and 
generalized limits (medium grey). In addition, we depict the outer limit of an e xtreme version 
of the generalized HZ (light grey) following the work by lMischna et al\ ( 120001 ). although this 
limit may be unrealistic based on subsequent studies. The orbit of HD 23079 b, a Jupiter- 
type giant planet, is depicted by a thick solid line. 
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Fig. 2. — Orbital stability simulations with HD 23079 b initially placed at periastron with 
= 1.503 AU, Cp = 0.071 and the Earth-mass planets placed at four different starting angles, 
which are: 45° (top left), 90° (top right), 270° (bottom left), and 315° (bottom right). Using 
a rotating coordinate system, HD 23079 b moves along the thin line. The Earth-mass Trojan 
planets, which give rise to the "banana-shaped" area at L4 or L5, remain within the HZ for 
at least 10^ years. 
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Fig. 3. — Orbital stability simulations with HD 23079 b initially placed at apastron with 
= 1.503 AU, Cp = 0.071 and the Earth-mass planets placed at four different starting angles, 
which are: 45° (top left), 90° (top right), 270° (bottom left), and 315° (bottom right). Using 
a rotating coordinate system, HD 23079 b moves along the thin line. The Earth-mass Trojan 
planets, which give rise to the "banana-shaped" area at L4 or L5, remain within the HZ for 
at least 10^ years. 
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Fig. 4. — Orbital stability simulations with HD 23079 b initially placed at apastron with 
= 1.503 AU, Cp = 0.102 and the Earth-mass planets placed at four different starting angles, 
which are: 45° (top left), 90° (top right), 270° (bottom left), and 315° (bottom right). Using 
a rotating coordinate system, HD 23079 b moves along the thin line. The Earth-mass Trojan 
planets, which give rise to the "banana-shaped" area at L4 or L5, remain within the HZ for 
at least 10^ years. 
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Fig. 5. — Orbital stability simulations with HD 23079 b initially placed at periastron for 
Earth-mass planets placed at starting angles of 180°. The model simulations differ regarding 
the selected values for the orbital parameters Up and Cp of the giant planet HD 23079 b 
(see also Table 2 for further information). The respective value pairs (a^, e^) are: (1.503, 
0.102), (1.596, 0.071), (1.596, 0.102), (1.596, 0.133), and (1.689, 0.102) for panel a, b, c, d, 
and e, respectively, with ap in AU. The grey domains (only visible in panel a and d) depict 
HD 23079's stellar HZ; see Table 2 for the times of first exit of the planet from the HZ. Note 
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Fig. 6. — Orbital stability simulations with HD 23079 b initially placed at periastron with 
ttp = 1.596 AU, Cp — 0.133 and the Earth-mass planets placed at a starting angle of 0°. 
The Earth-mass planet remains within the HZ for at least 10^ years. However, during that 
time it was captured by the giant planet and thus became a natural satellite (moon) of that 
planet, resulting in the small black area. Also note the absence of the "banana-shaped" area 
at L4 or L5. 
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Fig. 7. — Using the giant planet as reference and origin of the coordinate system, we measure 
the angle 9p of the captured terrestrial planet, which thus became a moon, in a sidereal frame 
(with = corresponding to the 3 o'clock position). With a uniform data sampling rate, 
the moon will be more likely recorded at or near apogee. Since the orbit of the moon is 
highly eccentric, there are many more points when the moon is near apogee compared to 
when it is near perigee. 
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Fig. 8. — Histogram displaying the time-dependent distance between the moon and the giant 
planet. Top: Simulation with an elapsed time of 1000 years and a data sampling period of 
0.01 years. Bottom: Simulation with an elapsed time of 1 x 10^ years and a data samphng 
period of 100 years. The comparison between both figures shows that there is little, if any, 
evolution of the system over the total time of simulation. 



- 22 - 



Table 1. Stellar and Planetary Parameters 



Parameter 


Value 


Reference 


Spectral Type 


F9.5 V 


Grav et al. f2006) 


RA 


3 39 43.0952 


ESA (1997)"''' 


DEC 


- 52 54 57.017 


ESA (1997)"''' 


Teff (K) 


6030 ± 52 


Ribas et al. f2003) 




1.106 ± 0.022 


Ribas et al. f2003) 


M{Mq) 


1.10 ± 0.15 


c 


Mv 


4.42 ± 0.05 


ESA (1997)"'''''^ 


Mboi 


4.25 ± 0.05 


ESA (1997)"'*'''^ 


Distance (pc) 


34.60 ± 0.67 


ESA (1997)"'''''^ 


Mpsini (Mj) 


2.45 ± 0.21 


Butler et al. f2006) 


P (days) 


730.6 ± 5.7 


Butler et al. f2006l 


S (AU) 


1.596 ± 0.093 


Butler et al r2006l 




0.102 ± 0.031 


Butler et al. f2006) 


Note. — "data 


from SIMBAD, 


see 



|http : / /simbad . u-strasbg . f r 



''adopted from the Hipparcos catalogue 

'^based on spectral type 

'^based on parallax 28.90 ±0.56 mas 
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Table 2. Time of First Exit (or Ejection) of the Earth-Mass Planet from the HZ (in years) 



ap 




0° 


45° 


90° 


135° 


180° 


225° 


270° 


315° 


(AU) 




















1.503 


0.071 


3.26E-01 U 






1.67E+02 L 


1.05E+02 L 


5.38E+01 L 






1.503 


0.102 


1.84E+00 C 


2.30E+05 U 


7.09E+05 U 


8.51E+01 U 


4.84E+01 L 


4.27E+01 U 


8.78E+05 U 


2.59E+05 U 


1.503 


0.133 


1.07E+00 L 


3.46E+02 L 


6.13E+02 L 


9.72E+01 U 


6.04E+01 U 


4.30E+01 U 


4.78E+02 U 


1.85E+03 U 


1.596 


0.071 


1.30E-03 C 


2.33E+02 U 


1.54E+01 U 


9.48E+00 U 


7.93E+00 U 


5.59E+00 U 


3.89E+00 U 


1.95E+00 U 


1.596 


0.102 


4.83E-01 U 


1.53E+01 U 


1.53E+01 U 


1.23E+01 U 


7.80E+00 U 


5.56E+00 U 


3.73E+00 U 


1.77E+00 U 


1.596 


0.133 




3.38E+01 L 


1.52E+01 U 


9.48E+00 L 


7.82E+00 U 


6.02E+00 U 


3.90E+00 U 


1.42E+00 C 


1.689 


0.071 


1.09E-01 C 


8.47E+00 U 


8.99E+00 U 


7.25E+00 L 


4.86E+00 U 


4.05E+00 U 


2.64E+00 U 


3.90E+01 L 


1.689 


0.102 


2.85E-02 C 


8.02E+00 C 


8.56E+00 L 


6.73E+00 U 


4.78E+00 U 


4.38E+00 L 


2.45E+00 U 


2.18E+00 L 


1.689 


0.133 


5.90E-03 C 


8.93E+00 U 


1.91E+01 U 


6.55E+00 L 


4.59E+00 U 


8.97E+00 U 


2.23E+00 U 


8.13E+00 U 



Note. — The total time of simulation is 10® yrs. The Jupiter-type planet started at the periastron position and the initial velocity of 
the Earth-mass planet was computed to begin a circular motion about the star. U means that the Earth-mass planet crosses the upper 
limit of the HZ given as 1.9662 AU; whereas L means that the Earth-mass planet crosses the lower limit of the HZ given as 0.9896 AU. C 

means that the Earth-mass planet has a close encounter with the giant planet, possibly resulting in a collision; therefore, the simulation was 
discontinued. If no data are given, the simulation lasted beyond 10® yrs without exiting the HZ. 



-24- 



Table 3. Time of First Exit (or Ejection) of the Earth-Mass Planet from the HZ (in years) 



Op 




0° 


45° 


90° 


135° 


180° 


225° 


270° 


315° 


(AU) 




















1.503 


0.071 


8.50E-01 U 






2.95E+01 U 


5.08E+01 L 


4.01E+01 U 






1.503 


0.102 


1.45E+00 U 






9.94E+01 L 


2.32E+01 L 


2.11E+01 U 






1.503 


0.133 


9.33E-01 L 


1.28E+03 L 


1.46E+03 U 


2.54E+01 L 


4.99E+01 U 


8.52E+00 U 


4.19E+03 U 


7.58E+02 U 


1.596 


0.071 


7.07E+00 L 


2.57E+02 U 


3.27E+01 L 


5.60E+01 C 


6.61E+00 U 


4.97E+00 U 


2.86E+01 U 


4.18E+01 U 


1.596 


0.102 


1.33E+01 U 


3.44E+01 U 


5.33E+01 U 


5.34E+01 U 


6.24E+00 C 


4.77E+00 U 


6.49E+01 U 


2.44E+01 U 


1.596 


0.133 


7.07E+00 U 


2.05E+01 U 


7.08E+01 U 


3.23E+01 U 


1.73E+01 U 


5.07E+00 U 


2.26E+01 U 


2.08E+01 U 


1.689 


0.071 


1.89E+01 C 


3.61E+01 U 


7.76E+00 U 


6.45E+00 U 


5.17E+00 U 


3.61E+00 U 


3.61E+00 U 


1.38E+00 U 


1.689 


0.102 


7.48E+00 L 


l.OOE+01 U 


7.59E+00 U 


7.32E+00 U 


5.52E+00 L 


3.16E+00 U 


3.32E+00 ET 


E2.SE+00 U 


l.()8<) 




1,87E+()1 i: 


i).()iE+0() r 


7.11E-()0 r 


7.02E-()0 E 


5.5:-!E+()l E 


:E27E+()0 E 


:E ilE+(X) E 


EO<sE+()0 U 



Note. — The total time of simulation is 10*" yrs. The Jupitci-typc planet started at the apastron position and the initial velocity of the 
Earth-mass planet was computed to begin a circular motion about the star. For further information see notes of Table 2. 



